Astronomy <fe Astrophysics manuscript no. SDSS-VLT 


© ESO 2008 


March 14, 2008 





Post common envelope binaries from SDSS. II: 
Identification of 9 close binaries with VLT/F0RS2* 

M.R. Schreiber^, B.T. Gansicke^, J. Southworth^ , A.D. Schwope'^, D. Koester^ 

^ Departamento de Fi'sica y Astronomia, Facultad de Ciencias, Universidad de Valparaiso, Valparaiso, Chile 
^ Department of Physics, University of Warwick, Coventry CV4 9BU, UK 
^ Astrophysikalisches Institut Potsdam, An der Sternwarte 16, D-14482 Potsdam, Germany 
* Institut fiir Theoretische Physik und Astrophysik, University of Kiel, 24098 Kiel, Germany 
e-mail: Matthias . SchreiberSuv . cl 

Received / Accepted 

Abstract 

Context. Post common envelope binaries (PCEBs) consisting of a white dwarf and a main sequence star are ideal 
systems to calibrate current theories of angular momentum loss in close compact binary stars. The potential held by 
PCEBs for further development of close binary evolution could so far not be exploited due to the small number of 
known systems and inhomogeneity of the sample. The Sloan Digital Sky Survey is changing this scene dramatically, 
as it is very efficient in identifying white dwarf/main sequence (WDMS) binaries, including both wide systems whose 
stellar components evolve like single stars and - more interesting in the context of close binary evolution - PCEBs. 
Aims. We pursue a large-scale follow-up survey to identify and characterise the PCEBs among the WDMS binaries that 
have been found with SDSS. We use a two-step strategy with the identification of PCEBs among WDMS binaries in the 
first phase and orbital period determinations in the second phase. Here we present first results of our ESO-VLT/FORS2 
pilot-study that has the target of identifying the PCEBs among the fainter ((;>18.5) SDSS WDMS binaries. 
Methods. From published SDSS catalogues we selected 26 WDMS binaries to be observed with ESO-VLT/FORS2 in 
service mode. The design of the ESO-VLT/FORS2 observations was to get two spectra per object separated by at least 
one night. We used the Nal AA 8183.27,8194.81 doublet to measure radial velocity variations of our targets and per- 
formed additional follow-up spectroscopy using Magellan-Clay/LDSS3 of two systems showing significant radial velocity 
variations. Using a spectral decomposition/fitting technique we determined the white dwarf effective temperatures and 
surface gravities, masses, and secondary star spectral types for all WDMS binaries in our sample. 

Results. Among the 26 observed WDMS binaries we find 9 strong PCEB candidates showing clear (> Scr) radial velocity 
variations, and we estimate the fraction of PCEBs among SDSS WDMS binaries to be 35± 12%. We find indications 
for a dependence of the relative number of PCEBs among SDSS WDMS binaries on the spectral type of the secondary 
star. These results are subject to small number statistics and need to be confirmed by additional observations. The 
orbital periods of two PCEB candidates, SDSS J1047+0523 and SDSS J1414-0132, we measured to be 9.17hrs and 
17.48 hrs respectively. 

Conclusions. This pilot study demonstrates that our survey is highly efficient in identifying PCEBs among the SDSS 
WDMS binaries, and will indeed provide the observational parameters that are necessary to constrain the theoretical 
models of close binary evolution. 
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1. Introduction 

The majority of all stars are born in binary or multiple 
star systems and a significant fraction of these will inter- 
act at some point during their lives. Once the more mas- 
sive stars leaves the main sequence and depending on the 
initial conditions, dynamically unstable mass transfer or a 
tidal instability may force the system to enter th e common 
envelope (CE) phase (see iTaam fc Rickeii l2006l for more 
details). In this phase a gaseous envelope around the entire 
binary forms and friction within the envelope is significantly 
reducing the binary separation. Since Paczviiski (,1976.) it 
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* Based on observations collected at the European Southern 
Observatory, Paranal, Chile under programme ID 078.D-0719. 
This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile. 



is generally believed that virtually all close compact bi- 
nary systems ranging from X-ray binaries to double white 
dwarf binaries or cataclysmic variables, to name a few, have 
formed through CE evolution. IWillems fc Kolbl (|2004l ) per- 
formed comprehensive binary population synthesis stud- 
ies for white dwarf/main sequence binaries (WDMS), and 
find that dynamically stable mass transfer may produce 
some additional short orbital period (Porb~ lOdays) sys- 
tems primarily consisting of low-mass He white dwarfs 
(0.1 < Mwd < O.4M0) and relatively massive secondary 
stars (up to 6M0). In the vast majority 75 — 99%) of 
cases, however, a CE phase is needed to produce short or- 
bital period (Porb^ 400 days) WDMS. 

From theoretical simulations we learned that the CE 
phase is probably very short, i.e. <10'^yrs, that the spi- 
ral in starts rapidly after the onset of the CE phase, and 
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that the expect ed form of post-CE planetary nebular are 
bipolar (see e.g. 'Morris'. '198lVTaam fc Bod enheimeH . ll989t 
[Taani & Sandquist, 2000; Wcbbink, _2007, and references 
theirin). Hydrodynamical studie s have been carried for a 
number of parameters (e.g. ITaam fc Rickeil . l2006l . for a 
recent review), however, these expensive simulations can- 
not be used in population synthesis models. Instead simple 
equations relating the total energy of the binary before an d 
after the CE phase are applied (e.g lWillems fc Kolbl . [200i) . 
Usually, the CE pha se is s i mply approxima ted by a para m- 
eterised energy (Wcbbinkl. [198 4: Willcms fc Kolhll2004D or 
angular momentum equation (fNelemans et al.l . 12000'). Both 
descriptions differ significantly in the predicted outcome of 
the CE phase, i.e. the energy equation predicts a much 
stronger shrinkage than the angular momentum equation 
(jNelemans fc Tout) . l2005l ). In addition, in both prescrip- 
tions the efficiency to "use" the orbital energy (or angular 
momentum) to expel the envelope is uncertain. Hence, the 
CE phase is probably the least understood part of close bi- 
nary evolution and currently intensivel y debated (fo r recen t 
discussion s of th e CE pha se see Nclc mans fc ToutI (j2005l) : 
iWebbinkI (|2007l) : iBeer et all (|2007D ). 

Once the envelope is expelled, the evolution of the 
PCEB is mainly driven by angular momentum loss 
(AML) due to magnetic braking. Unfortunately, the two 
currently favoured prescriptions for m agnetic braking 
(jVerbunt fc Zwaanl 119811: lAndronov eFall . |2003) , differ by 
up to two orders of magnitude in the predicted angular 
momentum loss rate. Even worse, it is not clear whether 
magnetic braking is present in all PCEBs, or whether it ter- 
minates in systems where the secondary star is fully convec- 
tive. In order to explain the orbital period gap observed in 
the period distribution of CVs, one in fact needs to assume 
the latter, known as "disrupte d magnetic braking", (e.g. 
iKingl 11988: How ell et"aLl . l200lD . In contrast to this, obser- 
vations of single low mass st ars do not show any evid ence 
for such a discontinuity (e.g. iPinsonneault et al.l . [2002[ ). 

Significant progress in our understanding of the CE 
phase and AML in PCEBs certainly requires both contin- 
uous theoretical efforts and innovative observational input. 
The potential that a large sample of white dwarf plus main 
sequence PCE Bs could provide in this con text has been 
first realised bv lSchreiber fc Gansickd ([2003') . They showed 
that it is possible to determine the orbital period the PCEB 
had when it left the CE phase and the white dwarf temper- 
ature as well as the orbital period it will have when entering 
the next phase of int eraction (becoming a cataclysmic vari- 
able). More recentlv JPolitano fc We ilcr (2006') showed that 
the disrupted magnetic braking hypothesis can be tested 
by knowing the relative number of PCEBs among white 
dwarf/main sequence (WDMS) binaries as a f u nction of 
secondary spectral type. Finallv [Schreiber et al.l ([2007t ) re- 
cently summarised that a representative sample of PCEBs 
could probably constrain the CE efficiency, the description 
of the CE phase, and the strength of magnetic braking. 

Unfortunately, only a very limited number of PCEBs 
with known effective temperature and orbital pe- 
riod is currently availabl e. The PCEB sample that 
[Schreiber fc Gansickd ([2003D analysed consisted of only 30 
well-studied short-period WDMS k nown in 2002. Since 
then. 
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of well studied PCEBs is close to 50. The number of sys- 
tems required to constr ain close b inary evolution is cer- 
tainly much higher. Sch reiber et a l. (2007) estimate that 
about 200 PCEBs with known white dwarf temperatures, 
secondary star parameters, and orbital periods are needed. 

The Sloan Digital Sky Survey (SDSS) is currently open- 
ing a new era of observational population studies for a wide 
variety of astronomical objects, including WDMS binaries. 
As by-product of the quasar search within SDSS, more than 
1000 WDMS bina ries have been id e ntified already within 
Data Release 5 ([Eisenstein et al] . [20061 [Silvestri et al.l 



|2007[: [Southwort h et al.[ . [20071 ). Though very large, this 



set of WDMS binaries is still subject to selection effects, 
as it favours systems with bluish colours, i.e. hot white 
dwarfs. To compensate for this, we have initiated a dedi- 
cated WDMS identification programme as part of SDSS II- 
SEGUE, which targets sy stems with a strong co ntribution 
of the companion star (see [Schreiber et al] . [20071 ) . In paral- 
lel, we have started the identification of PCEBs among the 
WDMS binaries identified by SD SS by means of a radial 
velocity survey. Rcbassa-Manscrgas et al] ([20071 henceforth 
Paper I) presented a first set of 18 PCEB candidates from 
radial velocity variations detected among multiple SDSS 
spectra. In this paper we present VLT/F0RS2 observations 
of 26 faint {g>18.5) WDMS binaries which led to the iden- 
tification of 9 new PCEBs through radial velocity varia- 
tions. The orbital periods of two of these PCEBs were de- 
termined from additional Magellan-Clay/LDSS3 observa- 
tions. The structure of the paper is as follows. In Section|5| 
we describe the selection of the 26 VLT targets and the 
determination of their stellar parameters. Details on the 
observations and data reduction of the VLT and Magellan- 
Clay runs are given in Sect.O After describing our method 
of measuring radial velocities and determining orbital pe- 
riods (Sect.dj) we discuss our results in the context of our 
large-scale PCEB survey. 

2. Target Selection and Stellar parameters 

For this project, we selected 26 WDMS binaries with i < 
19.5 and visible from Paranal from the SDSS Data release 



5 (DR5, [Adelman-McCarthy et al] [2007[). using the list s 



compiled bv Eisenstei n et 



tny et aij izuu/^ll. usmg tne nst s 
all (|2006D : [silvestri et all ((20071) : 



(2006), [Aungweroiwit et all ([20071 ). and iTappert et al 



2007f ) added additional systems and the current number 



Southwo rth et al, (2 0071 ) . Given the recently sugges t ed tes t 
of disrupted magnetic braking ([Politano fc Weile^ . [20061 ) . 
we focussed this selection on WDMS binaries with interme- 
diate M-type (M2-M5) companions. Tabled] Hsts the SDSS 
object names and the i magnitudes of our targets. 

We have determined the stellar parameters of these 26 
WDMS binaries following the method described in Paper I, 
with the only difference of now using spectra retrieved 
from Data Release 6, which have undergone a complete 
re-red uction that substantially improves the flux calibra- 
tion ([Adelman-McCarthv fc for the SDSS Coffaboratioi] , 
[2007[ ). For the sake of brevity, we provide here only a 
brief summary of the adopted procedures, for full details 
the reader is referred to Paper I. We first carry out a 
two-component fit to the SDSS WDMS spectra using 
libraries of observed white dwarf and M-dwarf template 
spectra, normalised to stellar surface fluxes. These tem- 
plate libraries are also compiled from SDSSDR6. This 
stage provides us with a spectral type for the companion 
star and a flux scaling factor for the M-dwarf template 
which, together with a spectral type-radius relation, yields 
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Figure 1. Two-component fits to the 2 PCEBs SDSS J1047+0523 and SDSS J1414-0132, for which we determined 
orbital periods from Clay/LDSS3 spectroscopy. Top panels: the WDMS spectra (black lines) along with the white dwarf 
and M-dwarf templates (dotted lines). The MJD~PLT-FIB identifiers of the two SDSS spectra are also given. Bottom 
panels: residuals from the fit. 
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Figure 2. Spectral model fits to the white dwarf components of SDSS J1047+0523 and SDSS J1414-0132, obtained by 
subtracting the best-fit M-dwarf template from their SDSS spectra. Top left panels: normalised H/3 to He line profiles 
(top to bottom, gray lines) along with the best-fit white dwarf model (black lines). Top right panels: 1, 2, and 3a 
contour plots in the T^g — logg plane. The black contours refer to the best line profile fit, the red contours to the fit of 
the whole spectrum. The dashed line indicates the location of maximum H/3 equivalent width in the T^g — logg plane, 
dividing it into "hot" and "cold" solutions. The best-fit parameters of the "hot" and "cold" normalised line profile 
solutions and of the fit to the 3850 — 7150 A range are indicated by black and red dots, respectively. Bottom panels: the 
white dwarf spectrum and associated flux errors (gray lines) along with the best-fit white dwarf model (black lines) to 
the 3850 7150 A wavelength range (top) and the residuals of the fit (gray line, bottom). 



a distance estimate. In a second step, we subtract the best- 
fit M-dwarf template from the WDMS spectra, and fit the 
residual white dwarf spectra with a standard method 
for variable white dwarf effective temperature (Tgff) and 
surface gravity (log 5), using a grid of pure-hydroge n (DA) 
model spectra computed with the code of .Koester et al.l 



(l2005l lP1. Two flavours of white dwarf fits were carried out 
for each object: (1) a fit to the normalised line profiles. 



We limit the spectral fitting to DA white dwarfs, but provide 
a spectral type of the companion for the two helium atmosphere 
(DB) WDMS binary in our sample. 
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Figures. Distances of the 24 WDMS binaries in Table[T] 
that contain a DA white dwarf. The spectral modelling 
(Sect. [3]) provides two independent distance estimates for 
the individual stellar components. Nine systems have dgcc > 
dwd by more than 1 a. Adopting secondary star spectral 
types later by one class brings most of them onto the 
dscc — f^wd relation. An earlier appearance of the spectral 
type determined from an optical spectrum may be related 
to stellar activity. See Paper I for additional discussion. 



Due to the degeneracy in Balmer line equivalent width as 
a function of temperature two acceptable fits are usually 
found, a "hot" and a "cold" solution. (2) a fit to the 
continuum plus lines over the range 3850 — 7150 A. The 
parameters from this fit are subject to substantial system- 
atic uncertainties due to inaccuracies in the spectral slope, 
and we only use them to break the degeneracy between 
the "hot" and "cold" line profile solutions. Examples of 
the spectral decomposition and of the white dwarf fits are 
shown in Fig. [1] and Fig. (5] 

We then used an updated version of the tables 
in [Bergeron et al.l ()1995f n. which inc orpora te the struc- 
ture and evo l ution calculations of IWoodI (|1995f ) and 
I Fontaine et al.l (j200lh , to compute white dwarf masses and 
radii from Tcff and logg. Finally, the white dwarf radii and 
flux scaling factors of the white dwarf models yield a second 
distance estimate. Statistical errors on the different prop- 
erties are propagated through this process, but should be 
considered with some caution as additional systematic un- 
certainties are not easy to quantify. The resulting stellar 
parameters are listed in Table[T] Inspection of Fig. [3] (left 
panel) shows that in about one third of the systems the dis- 
tance estimate based on the fit of the secondary star exceeds 
that based on the white dwarf fit by a significant (> 1 cr) 
factor. A similar ratio of systems with dscc > c'wd has been 
found in the much larger sample analysed in Paper I, which 
also discusses in detail possible explanations for this find- 
ing. Our favoured hypothesis from Paper I is that stellar 
activity causes large scatter in the spectral type-radius re- 
lation of low-mass stars (in fact, larger than the scatter in 
their mass-radiu s relation that is subject to thor ough in- 
vestigations, e.g. lRibasll2006t lL6pez-Moralesll2007^ . with a 
tendency for spectral types too early for their radius. In 
the right-hand panel of Fig. [31 we show that adopting sec- 
ondary stars of one spectral class later than in Table [T] 
brings the nine systems which had c?soc > c'wd on, or close 
to, the dscc = c?wd relation. 

A final note concerns SDSS J2324-0931, which con- 
tains a DB white dwarf. The DR3 SDSS spectrum of 

^ http://www.astro.umontreaI.ca/ bergeron/CoolingModels/ 



SDSS J2324-0931 was analysed by| van den Besselaar et al.l 
(|2005f ) who adopted a method similar to our spectral de- 
composition, i.e. fitting the WDMS spectrum with the com- 
bination of observed white dwarf and M-dwarf templates. 
Different to our approach, they fixed the system distance 
to be the same for both stars. Their analysis resulted in 
Toff = 36815, Sp(2) = M3V, and d = 911 pc. Our spectral 
decomposition results in Sp(2) = M4V and d = 633±188pc. 
Forcing the spectral type to M3V, we obtain a distance 
of 879 ± 140 pc. Given that the typical uncertainty in the 
spectral type is one spectral class, the two results are hence 
consistent. 

3. Observations 

3.1. VLT/F0RS2 

Intermediate resolution spectroscopy of 26 WDMS bina- 
ries were obtained between August 16, 2007 and March 
26, 2007 with F0RS2 on the VLT/UTl (Antu). The ex- 
posure time was 900 sec for all objects. The observations 
were carried out using the 1028z grism and a 1" slit, re- 
sulting in a spectral coverage of 7830-9570 A. The images 
were processed using the STARLINK packages FIGA RO and 
KAPPA, and the spectra were op timally extracted (jHornel 
1986) using the PAMELA package (lMarshl. ll989t ). From mea- 
surements of the FWHMs of sky lines we find our observa- 
tions have a resolution of approximately 2.5 pixels (2.2 A) at 
8200 A. Wavelength calibration of the extracted spectra was 
done using only sky emission lines. The wavelengths o f good 
sky lines we re obtained from the atlas of Osterbrock et al.l 
(|l996l . ll997i l. 37 sky lines were included, and fitted with a 
fifth-order polynomial. The rms scatter around the fit was 
0.11 A, s o the statistical unce rtainty in the wavelength scale 
is 0.04 A (jMarsh et al.l .'l994'). Finally, the spectra were flux- 
calibrated and compensated for telluric absorption using a 
spectrum of the standard star FeigellO. 

The Phase II design of these observations was to obtain 
two spectra per object, separated by at least one night, in 
order to probe for radial velocity variations. Due to poor at- 
mospheric conditions, additional repeat observations were 
obtained for several targets, while two objects only got one 
F0RS2 spectrum each. 

3.2. Magellan-Clay/LDSS3 

Follow-up spectroscopy of two PCEBs identified from the 
F0RS2 data, SDSS J1G47+0523 and SDSS J1414-G123, 
were obtained over the period 2007 May 17-20 using the 
Magellan Clay telescope equipped with the LDSS3 imaging 
spectrograph. Seeing and transparency were highly vari- 
able. We adjusted the exposure times accordingly, ranging 
from 500-900sec for SDSS J1047+0523 and 750-1200sec 
for SDSSJ1414~0123. We used the VPH.Red grism and 
OG 590 order-sorting filter. The detector was an unbinned 
STA 4kx4k pixel CCD read out by two amplifiers. A 0.75" 
offset slit gave a wavelength coverage of 5800-9980 A at 
a reciprocal dispersion of 1.2Apx^^. We took flat-field 
images at each sky position whilst observing as this was 
needed to remove the effects of fringing at redder wave- 
lengths. The dispersed images were processed using the 
same software as for the VLT observations. Measurements 
of the FWHMs of the sky lines indicate that our obser- 
vational setup gave a resolution of approximately 4 pixels 



Schreiber et al.: Identifying PCEBs with VLT/F0RS2 



5 



Table 1. Stellar parameter of the 26 WDMS binaries observed with F0RS2, as determined from spectral modelling 
(Sect. 121). Upper limits on the orbital periods for the nine PCEBs identified by our VLT/F0RS2 observations are calculated 
as described in Sect. 15. II 



System (SDSS J) MJD-plate-fibre 



r[Ki 



log 3 Afwd[MQ] dwd[pc] Sp dscc[pc] 



notes 



001726.63-002451.1 



022835 
032038 
081959. 
083354. 
085542. 
092451 
093904. 
095108. 
100953 
104358. 
104738. 
110151. 
112623. 
122208. 
125250 
141451 
142541. 
150657, 
204541, 
205059, 
212320. 
215614 
220338 
231221. 
232438 



.92-074032.3 
.72-063822.9 
20+060424.2 
84+070240.1 
49+044717.7 
.63-001736.4 
03+051114.8 
74+025507.5 
.69-002853.4 
59+060320.9 
24+052320.3 
34+122241.5 
87+010856.8 
48+093406.2 
.03-020608.1 
.73-013242.7 
32+023047.4 
58-012021.7 
.90-050925.7 
.37-000254.3 
74+002455.5 
.57-000237.4 
.61-001750.7 
59+010127.0 
.31-093106.4 



52559-1118-280 
52518-0687-153 
51908-0454-534 
51924-0460-432 

52962- 1296-138 

52963- 1297-511 
52670-1190-512 
52000-0474-175 
52710-0993-296 

51908- 0481-419 

51909- 0270-181 
52643-1000-055 
52339-0578-609 
53119-1603-492 
51614-0281-402 
52672-1230-73 
51694-0338-343 
52400-0917-329 
51999-0535-348 
52426-0922-536 
52145-0635-439 
52466-0982-567 
52523-0987-484 
52078-0371-025 
52173-0372-034 
51811-0381-570 
52203-0645-637 



18.3 
18.3 
18.7 
18.7 
19.2 
18.5 
19.1 
17.6 
19.2 
18.5 
19.0 
18.6 
18.7 
18.8 
17.5 
18.3 
18.5 
19.0 
18.9 
18.4 
18.7 
19.2 
19.2 
18.5 
18.7 
17.3 
18.4 



12828+ 2543 
13588+ 1803 
13588+ 948 
11045+ 479 
15782+ 1422 
16336+ 730 
10073+ 336 
24163+ 1359 
12681+ 1406 
13432+ 1007 
21535+ 1927 
11565+ 536 
12392+ 1847 
20331+ 817 
16717+ 781 
30419+ 2053 
22037+ 2333 
14065+ 1452 
13127+ 1014 
15782+ 1005 
26801+ 1432 
18330+ 1169 
12250+ 1794 
7910+ 926 
32972+ 1703 



8.00+0.46 
8.11+0.39 
7.76+0.17 
8.49+0.33 
7.67+0.34 
8.03+0.16 
8.46+0.43 
7.82+0.20 
7.69+0.42 
7.80+0.22 
8.35+0.32 
8.39+0.24 
7.54+0.42 
7.83+0.15 
7.91+0.18 
7.34+0.49 
8.29+0.38 
8.31+0.27 
7.91+0.22 
7.72+0.24 
7.78+0.22 
8.27+0.24 
7.47+0.42 
9.50+0.74 
7.24+0.38 



0.61+0.29 
0.68+0.25 
0.48+0.09 
0.91+0.19 
0.45+0.17 
0.63+0.10 
0.89+0.25 
0.54+0.11 
0.44+0.23 
0.50+0.12 
0.84+0.19 
0.85+0.15 
0.38+0.20 
0.53+0.08 
0.57+0.10 
0.38+0.18 
0.80+0.23 
0.80+0.17 
0.56+0.13 
0.47+0.13 
0.53+0.11 
0.78+0.15 
0.35+0.19 
1.37+0.21 
0.37+0.11 



423+119 
423+107 
423+ 43 
247+ 61 
764+150 
424+ 45 
247+ 77 
688+ 89 
591+139 
425+ 57 
609+141 
220+ 39 
474+111 
534+ 50 
505+ 57 
1891+599 
641+171 
402+ 77 
362+ 48 
584+ 85 
892+127 
539+ 93 
496+119 
86+ 76 
1772+477 



577+171 
529+157 
311+140 
461+232 
860+255 
666+198 
589+297 
588+ 94 
919+273 
726+215 
867+257 
256+115 
508+256 
1071+171 
733+152 
1543+166 
1282+265 
552+278 
390+176 
629+187 
1007+160 
998+296 
606+305 
346+174 
1598+379 
455+ 72 
633+188 



1, Porb < 143.1 d 



-Porh 



< 14.7d 

< 61. 9d 



1, 2, Port < 0.8 d 



1, 3, 6, Porb < 3.4 d 

1, Porb < 9.1 d 
1, Porb < 49.1 d 

1, Porb < 0.4 d 

1, 4, 6, Porb < 44.0 d 

5 
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(1) PCEB candidate on the base of radial velocity variations detected among the available FORS2 and SDSS spectra; (2) 
Porb ~ 9.17 h determined from Magellan-Clay/LDSS3 spectroscopy; (3) Porb = 17.48 h determined from Magellan-Clay/LDSS3 
spectroscopy; (4) blue component is faint, white dwarf parameter uncertain; (5) DB white dwarf; (6) poor quality spectrum. 



(4.8 A) at 8200 A. Wavelength calibration was done using 
only sky lines, to increase observational efficiency. A total 
of 36 sky lines were fitted by a fifth-order polynomial. The 
rms scatter around the fit was 0.25 A, so the statistical un- 
certainty is 0.09 A. From this work and other studies we 
have found that reliable wavel ength calibrations can b e ob- 
tained from sky hues (see also lSouthworth et al.ll200^ . but 
that the scatter of the fit around the wavelength solution is 
quite a bit larger than for arc-lamp emission lines. In this 
case the statistical uncertainty is equivalent to 0.09 A so 
the increased scatter is relatively unimportant. 

Flux calibration and telluric line removal was performed 
using spectra of the standard star LTT3218 obtained dur- 
ing the same observing run. 

4. Radial velocity variations and orbital periods 

We have used the strong Nai AA 8183.27,8194.81 absorption 
doublet present in the spectra of all our targets to measure 
the radial velocities of the main-sequence companion stars. 
For this purpose, we fitted the wavelength of the Nai dou- 
blet using the sum of a linear slope for the continuum plus 
two Gaussians for the absorption lines. We fixed the separa- 
tion of the two Gaussians to 11.54 A corresponding to the 
laboratory value. Free parameters in the fit were the con- 
tinuum flux level and slope, the wavelength of the doublet, 
the full- width at half maximum (FWHM) of the lines (both 
FWHM were fixed to vary together), and the amplitudes 
of the individual doublet components. The fits were carried 




500 
Velocity [km s"' 



Figure 5. Examples of the double-Gaussian fits 
to the Na I AA 8183.27,8194.81 line profiles in the 
SDSS, VLT/F0RS2, and Clay/LDSS3 spectra of 
SDSS J1047+0523. 
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Figure 4. Radial velocities (Table[3] of our sample 26 WDMS binaries (Table[T]), measured from the 
Na I AA 8183.27,8194.81 doublets observed in their SDSS and F0RS2 spectra. No SDSS radial velocities could 
be determined for for SDSS J1252-0206, SDSS J1414-0132, and SDSS J2050-0002. See Fig.lH for examples of the 
Nai AA 8183.27,8194.81 fits. The systems with black symbols display radial velocity variations at a 3cr level. 
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Figure 6. Scargle periodograms calculated from the 
Na I AA 8183.27,8194.81 radial velocities measured in the 
LDSS3 spectra of SDSS J1414-0123 and SDSS J1047+0523. 

out using the fit/table command in MIDAS. Examples of 
these fits are shown in Fig.[5l We added the uncertainty 
in the wavelength calibration in quadrature to the statisti- 
cal errors in the wavelength of the Nai resulting from the 
double-Gaussian fit. 

One system, SDSS J1222-h0934, has an MOV secondary 
star, which exhibits no significant Nai absorption. For this 



single object, we used the same method as above, but fitted 
three Gaussian lines of fixed separation and common width 
to the CaiiAA 8498.02,8542.09,8662.14 triplet. 

We combined the radial velocities obtained from the 
F0RS2 data with those measured from the DR6 SDSS 
spectra (Table[3]). Comparing the individual radial veloc- 
ities of each object, we define as a strong PCEB candi- 
date all systems which show a radial velocity variation at 
a > 3(7 level. Applying this criterion, we identify 9 PCEB 
candidates among our target sample of 26 WDMS binaries 
(Fig.a TablelB). 

The radial velocity variations of the companion stars 
in SDSSJ1414-0132 and SDSS J1047-h0523 were deter- 
mined in an analogous fashion as described above from 
the C lay/LDSS3 spectroscopy. We then calculated Scargle 
(|1982D periodo grams from the mean-subtracted radial ve- 
locity measurements to establish the periods of the two 
WDMS binaries (Fig.[S]). We carried out sine-fits to the 
radial velocities using the three highest aliases in the pe- 
riodograms as initial value for the frequency, and report 
the resulting periods, 7-velocities and radial velocity am- 
plitudes, as well as their errors and the reduced values, 
in TableU 

The lowest reduced values as well as the lowest rel- 
ative errors on the fit parameters are found for the periods 
corresponding to the highest peaks in the periodograms. We 
therefore conclude that the orbital periods of SDSS J1414- 
0132 and SDSS J1047-I-0523 are Porb = 17.48 ± 0.07 h 
and Porb = 9.17 ± 0.02 h. Folding the radial velocities 
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from the SDSS spectrum clearly identifies the system as a 
close binary. 



5. Discussion 

5.1. Orbital period limits 

As in lRebassa-Mansergas et all (|2007D we determine an up- 
per limit for the orbital period based on the measured ra- 
dial velocity variations (Table[31)- The orbital period as a 
function of the stellar masses, the binary inclination (i), 
and the radial velocity amplitude (i^scc) of the secondary 
follows from Kepler's 3rd law: 



(Afwd + Msoc)^ 



2ttG 



(1) 



0.5 



1 1.5 
Orbital Phase 



Figure 7. The NaiAA 8183.27,8194.81 radial veloci- 
ties of the companion stars in SDSS J1414-0132 and 
SDSS J1047-I-0523 folded over their respective orbital pe- 
riods of Porb = 17.48h and Porb = 9.17h. 

Table 2. Orbital periods, semi- amplitudes, and re- 
duced from sine-fits to the radial velocity data of 
SDSS J1047-h0523 and SDSS J1414-0132 for the strongest 
three aliases in the periodograms shown in Fig. [51 The best- 
fit values are set in bold font. 



System 



Porb [h] K2 [km s- 1 ] 7 [km s- 1 ] 



SDSSJ1047+0523 


15.11±0.12 


195±16 


20±11 




9.17±0.02 


187±7 


45±5 




6.64±0.04 


181±22 


66±17 


SDSSJ1414-0132 


51.01±0.53 


371±32 


44±12 




17.48±0.07 


168±4 


34±3 




10.14±0.09 


152±13 


46±10 



over the neighbouring aliases results in substantially dis- 
torted non-sinusoidal periods. The LDSS3 radial veloci- 
ties of SDSS J1414-0132 and SDSS J1047-H0523 folded over 
their orbital periods are shown in Fig.[7l 

Below we provide some comments on individual objects. 



SDSS J0017-0024: Two SDSS spectra are available for this 
object, with radial velocities that are consistent at a 2a level 
()Rebassa-Mansergas et"all . l2007[ ). We obtained two F0RS2 
spectra, which agree within their la errors, confirming this 
system as a likely wide binary WDMS. 



SDSSJ1047-h0523: Only a single VLT/F0RS2 spectrum 
was obtained of this object. However, combining the ra- 
dial velocity from the F0RS2 spectrum with that measured 



Using the stellar parameters estimated in Sect. [2] from 
fitting the SDSS spectrum, assuming i = 90°, and as- 
suming the largest radial velocity difference measured 
for one particular system to be equal 2 Kscc , clearly 
gives an upper limit on Pprb- This method was used in 
iRebassa-Mansergas et al.l ()2007t ). Here we slightly improve 
this approach by using also the information provided by the 
times of observations. If the time delay between two spectra 
(At) is shorter than half the considered orbital period, the 
maximal radial velocity difference that could have been ob- 
served writes 2ifsoc sin(7rAi/2Porb) and we modify Eq. (1) 
accordingly. Upper limits on the orbital period for the 9 
PCEBs obtained with this method are given in Table[T] 

5.2. The fraction of PCEBs among SDSS WDMS binaries 

The first result of our PCEB survey will be the relative 
number of PCEBs among SDSS WDMS binaries. According 
to th e derived limits on the orbital periods and the masses 
g^en in Table 1 , the 9 close WDMS binaries identified in 
— tliia paper clearly experienced a phase of mass transfer dur- 
the post main sequence evolution of the progenitor of 
white dwarf. As mentioned in the Introduction, ac- 
^l^(grding to current WDMS binary population studies (with 
■jii^ny uncertainties involved), the majority of all WDMS 
l^^t formed through mass transfer interactions are ex- 
pBcted to be PCEBs and only a small fraction of close sys- 
tems may form through dynamically stable mass transfer. 
For the 9 PCEB candidates identified among our 26 WDMS 
systems dynamically stable mass transfer appears to be 
extremely unlikely as the white dwarf masses are rather 
high and the secondar y spectral types are rather late (see 
IWillems fc Kolbl . 12004 , for more details). We therefore state 
to have identified 9 PCEBs. 

In the first paper of this series 
(jRebassa-Mansergas et al.L 120071 ) multiple SDSS spec- 
troscopy has been used to identify 18 PCEB candidates 
among 101 WDMS binaries from the SDSS. Due to the low 
spectral resolution of SDSS spectra and the fact that in 
most of those 101 WDMS binaries only two SDSS spectra 
are available, the resulting PCEB fraction of ~ 15% 
certainly has to be considered as a lower limit on the true 
value. In this paper we identified 9 PCEBs among 26 SDSS 
WDMS binaries which gives a PCEB fraction of 35 ± 12%, 
which is - as expected - substantially higher than our 
previous lower limit. The increased PCEB detection rate 
is probably the result of both, the higher resolution of 
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VLT/F0RS2 and the fact that for most systems 3-4 
spectra are available. 

However, also in the case that more than 2 spectra are 
available exists a certain probability to sample similar or- 
bital phases and also the observations presented here have 
a finite spectral resolution. In the context of our project on 
PCEB studies, it is important to understand possible se- 
lection effects in the established PCEB sample. We there- 
fore developed a Monte-Carlo code to estimate the detec- 
tion probability among SDSS WDMS binaries as a func- 
tion of orbital period. For the purpose of this specific pa- 
per, we limit our analysis to the following assumptions: a 
mean radial velocity resolution of 6km/s, stellar masses 
of Mwd = 0.6 M© and Afgoc — 0.3 Mq, and uniformly dis- 
tributed inclinations. For an assumed orbital period we ran- 
domly select the inclination, orbital phase, and times of 
observations. The only constraint on the latter being that 
the time differences should exceed 12 h and that all observa- 
tions are performed in one semester. The randomly selected 
values are then used to calculate radial velocity differences. 
Repeating this exercise 10** times and counting the num- 
ber of cases with a Scr (alternatively 1 or 2 cr) radial ve- 
locity variation gives the probability of detecting a PCEB 
with the assumed orbital period. Fig. [5] shows the results 
we obtain assuming that 2 or 3 spectra have been taken 
in one semester. In both cases the detection probability 
of different significance levels of radial velocity variations 
is shown. Of course, our strategy leads to a certain bias 
towards shorter orbital periods (and higher inclinations) 
because strong radial velocity variations are more likely 
to be detected in these systems. However, with three ran- 
domly taken VLT/F0RS2 spectra our 3(t criterion is still 
~ 15 — 30% sensitive to PCEBs with periods in the range of 
^ 20 — 60 days. PCEBs with orbital periods longer than 70 
days may hide among the systems with weak (2cr) radial ve- 
locity variations, and will require additional spectroscopy 
to be unambiguously identified. In general. Fig. [5] clearly 
shows that the fraction of PCEBs among SDSS WDMS bi- 
naries of 35 ± 12% obtained here still represents a lower 
limit. However, as mentioned in the Introduction, the rel- 
ative number of PCEBs may be a function of the spectral 
type of the secondary. 

5.3. Dependence on the secondary star spectral type 



Recently, iPolitano fc Weileil (|2006D suggested a test for the 
disrupted magnetic braking scenario, i.e. the hypothesis 
that magnetic braking is not active in PCEBS which con- 
tain a fully convective companion star. Their calculations 
predict that the relative number of PCEBs among WDMS 
binaries should drastically increase for secondary star 
masses below M^cc = O.35M0. Our results are consistent 
with that prediction. Among 9 WDMS binaries with spec- 
tral type M3 or earlier we find only one strong PCEB can- 
didate while the PCEB identification rate among WDMS 
binaries containing later secondaries (M4-M6) is 8/17. 
According to the Poisson statistic, the relative numbers 
of PCEBs among SDSS WDMS binaries are 0.111 ± 0.111 
(for M0-M3) and 0.471 ± 0.166 (for M4-M6). Combining 
these values we obtain a significance of 1.8a for the hy- 
pothesis of having two different distributions. Interestingly, 
a potential difference between the secondary star spectral 
ty pes in close and wid e WDMS binaries was also suggested 
by iFarihi et al.l ()2005[ ). who found, subject to small num- 




Figure 8. Detection probability versus orbital period for 
2 and 3 spectra taken randomly in one semester but not 
in the same night. We here assumed Mwd = O.6M0 and 
Msec = O.3M0. The l — 3a labels refer to the detection level 
of radial velocity variations. Apparently, using our 3(T crite- 
rion to select strong PCEB candidates for further follow-up 
observations leads to final a sample that is biased toward 
shorter orbital periods. However, only 3 spectra are needed 
to identify 15 — 30% of the PCEBs with longer orbital peri- 
ods (Porb~ 20 days) and our strategy can therefore certainly 
be used to also constrain the orbital period distribution of 
long orbital period PCEBs. 



ber statistics, that the median secondary spectral type of 
close binaries is one spectral class later than that of the 
wide binaries. While keeping in mind the well known prob- 
lems resulting from low number statistics, these results may 
represent first indications for a discontinuity in the rela- 
ti ve number of PCEBs amo ng WDMS binaries as predicted 
bv lPolitano fc Weileil (|2006[ ). However, further observations 
are certainly needed to confirm or disprove the observed 
tendency. 

5.4. Predicted evolution of SDSS J1047+0523 and 
SDSS J1414-0132 



As shown by ISchreiber fc Gansickd (|2003l ) knowing the 
mass and effective temperature of the white dwarf, its or- 
bital period, and at least reasonable estimates for the mass 
and the radius of the secondary one can determine the age 
of the PCEB, the orbital period of the future semi-detached 
system (i.e. a cataclysmic variable) and the time the PCEB 
evolution takes for a given system. 

To that end, we here first follow [Schreiber fc Gansicke 
(.20031) by using the white dwarf cooling tracks of IWood 
()1995D to determine the white dwarf cooling ages of 
SDSS J1047-)-0523 and SDSS J1414-0132. Using the white 
dwarf parameters given in Table 1 and interpolating the 
cooling tracks we get icooi ~ 2.1 ± 1 x lO^yrs and tcooi 
3.9 ± 1.1 X 10**yrs for SDSS J1047-H0523 and SDSS J1414- 
0132 respectively. 

The determination of the mass and the radius of the 
secondary is much more challenging as the inclination of 
these two systems is unconstrained. We use the empir- 
ical spectral type-radius-mass (Sp-R-M) relation derived 
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Figure 9. The PCEB evolution of SDSS J1047+0523 and 
SDSS J1414-0132. Starting at their current orbital periods 
of Porb = 17.48h and Poib — 9.17h the systems evolve 
towards shorter orbital period due to angular momen- 
tum loss by gravitational radiation. The companion stars 
in SDSS J1047+0523 and SDSS J1414-0132 will fill their 
Roche-lobes in 1.7 ± 0.3 x 10^° yrs and 1.8 ± 0.2 x 10" yrs, 
respectively, and turn these PCEBS into cataclysmic vari- 
ables. This will happen at very similar orbital periods, 
0.090 ± 0.013 days for SDSS J1047-I-0523 and 0.085 ± 0.015 
for SDSS J1414-0132. 



6. Summary 

We have obtained ES0-VLT/F0RS2 spectroscopy of 26 
WDMS binaries found in SDSS. Clear radial velocity vari- 
ations led to the identification of 9 strong PCEB candi- 
dates, which gives a hit-rate of 35 ± 12%. We analysed 
the PCEB detection probability of our strategy and con- 
clude that (although a bias towards short orbital period 
systems is inevitable) our survey is sensitive to the pre- 
dicted long orbital period distribution of PCEBs. While 
we find 8 PCEB candidates among 17 systems with M4 or 
later secondary stars only one PCEB candidate could be 
identified among 9 WDMS binaries containing earlier sec- 
ondary stars. Although suffering from low number statis- 
tics, one may interpret this results as a first indication for 
a discontinuity in the dependence of the relative number of 
PCEBs among WDMS binaries on the spectral type of the 
secondary. 

We determined the orbital periods for two of the PCEBs 
identified here, SDSS J1047+0523 and SDSS J1414-0132, 
from Magellan-Clay/LDSS3 spectroscopy to be to be 9.17h 
and 17.48 h respectively. Both systems will enter the semi- 
detached CV state in or below the orbital period gap. This 
evolution will last longer than the Hubble time and both 
systems are therefore different than the progenitors of the 
current CV population. 

Acknowledgements. MRS acknowledges support from FONDECYT 
(grant 1061199), DIPUV (project 35), and the Center of Astrophysics 
in Valparaiso. 



in Paper I (Sect. 3.4) to estimate values of the mass and 
radius of the secondary. Apparently, the scatter around 
this empirical relation is quite large (Paper I, Fig. 7) and 
therefore the obtained values should be generally consid- 
ered as uncertain. However, in the case of SDSS J1414-0132 
and SDSS J1047+0523 the distance estimates based on the 
white dwarf (dwo) and those derived using the empirical 
Sp-R-M relation agree quite well which suggests that the 
obtained values are plausible (for a detailed discussion of 
the spectral type-mass-radius relation see Paper I). As both 
secondaries are of the spectral type M5 the empirical rela- 
tion gives Msec = O.2I6M0 and i?scc = 0.258^0. For our 
theoretical analysis we adopt rather broad ranges of pos- 
sible values to account for the described uncertainties, i.e. 
Msec = 0.20 - 0.23Mq and Psec = 0.2 - O.3i?0. 

According to the standard scenario of CV evolution, 
magnetic braking gets disrupted when the secondary stars 
becomes fully convective (e.g. IVerbunt fc Zwaaiil |1981[ ) 
and, hence, gravitational radiation is the only angular mo- 
mentum loss mechanism acting in systems with a low 
mass (Mscc^O.SMq) secondary star. The evolution of 
SDSS J1047-I-0523 and SDSS J1414-0132 according to this 
scenario and taking into account the uncertainties discussed 
above, is shown in Fig. [51 Both systems will become CVs in 
or slightly below the 2 — 3h orbital period gap and for both 
systems the evolution will take longer than a Hubble-time. 
Therefore, according to Schreiber & Gansickc (2003), they 
are not representative for progenitors of the current CV 
population. 
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Table 3. Radial velocities measured from the 
Na I AA 8183.27,8194.81 doublet in the ESO/VLT and 
SDSS spectra. The radial velocities obtained from SDSS 
data is indicated by the HJD set in italics. 



HJD RV [kms-i] HJD RV [kms-^] 



SDSS J0017-0024 

2453978.64799 -26.7±5.0 
2453982.67389 -18.3±5.5 
2452518.922000 -1.8±16.7 
2452559.785000 -36.3±15.5 

SDSS J0228-0740 
2453983.68162 4.0±8.0 
2453983.69278 9.9±7.5 
-16.9±21.0 

SDSS J0320-0638 
2454030.85227 240.0±5.9 
2454086.60870 204.1±7.0 
2454124.60190 212.1±6.0 
2451924.67646 172.9±24.8 

SDSS J0819+0604 
2454030.87927 26.4±5.7 
2454055.83644 39.2±6.0 
2452963.00286 19.8±17.1 

SDSS J0833+0702 
2454055.85156 83.5±5.5 
2454059.82430 121.6±5.9 
2452963.98436 194.5±18.3 

SDSS J0855-h0447 
2454059.83910 73.3±7.3 
2454060.84131 -11.7±7.0 
2452670.87671 15.8±21.1 

SDSS J0924-0017 
2454067.842758 15.4±4.5 
2454086.796712 14.7±5.2 

SDSS J0939+0511 
24.54086.81276 103.7±6.8 
2454087.84233 86.5±5.4 
2452710.81094 130.1±24.9 

SDSS J0951+0255 
2454086.82795 13.6±7.4 
2454099.81594 20.1±5.6 
2451908.91909 37.0±18.3 

SDSS J1009-0028 
2454086.84235 -2.6±9.0 
2454112.83168 -2.6±8.8 
2451909.91534 -48.0±23.7 

SDSS J1043+0603 
2454117.75282 46.2±7.0 
2454149.70110 48.0±7.2 
2452644.01881 37.7±22.9 

SDSS J1047+0523 
2454142.77188 -19.4±8.9 
2452339.82994 218.0±19.8 

SDSS JllOl+1222 
2454142.78686 -16.1±9.9 
2454155.63799 -22.0±8.8 
2454167.69909 -23.1±6.9 
2454169.67087 -25.3±6.0 
2454184.61366 -40.3±6.3 
2453119.74299 -56.8±19.8 



SDSS J1126+0108 

2454142.80234 -7.7±4.9 
2454155.62142 -12.1±4.7 
245I6I4.8O8I4 -1.1±15.4 

SDSS J1222+0934 
2454142.831916 -50.4±4.6 
2454156.875771 -48.3±4.5 

SDSS J1252-0206 
2454142.84647 -30.4±6.4 
2454169.68656 -26.7±5.6 
2454185.58421 -31.1±6.0 
2451694.70928 -33.3±24.0 

SDSS J1414-0132 
2454142.86192 154.2±7.3 
2454169.71281 -65.9±5.7 
2454180.89912 137.7±7.8 
35.2±24.9 

SDSS J1425+0230 
2453982.48632 102.6±8.4 
2453983.49473 95.3±19.7 
2451999.89365 101.1±26.3 

SDSS J1506 0120 
2453965.51719 41.8±5.3 
2453982.50834 -85.4±5.6 
2452426.79507 61.2±18.6 

SDSS J2045-0509 
2453978.52649 -48.4±5.8 
2453978.53830 -42.5±5.9 
2453981.61523 -21.6±5.2 
2452145.71287 12.1d=21.2 

SDSS J2050-0002 
2453998.625440 34.1±9.1 
2454001.522929 15.0±5.8 
2452466.81829 30.8±37.5 

SDSS J2123+0024 
2453978.55706 56.1±14.5 
2453981.63252 -91.6±7.8 
2452523.69331 198.2±35.6 

SDSS J2156-0002 
2453992.53760 -18.0±5.4 
2453998.64315 -98.9±6.4 
2454004.55526 21.6±5.1 
2452075.89855 -68.1±15.9 

SDSS J2203-0017 
2453978.57300 -48.7±6.4 
2453982.53714 -61.9±7.2 
2452173.63031 -20.9±29.5 

SDSS J2312+0101 
2453978.59006 14.7±4.7 
245I81 1.73493 18.3±16.2 

SDSS J2324-0931 
2453963.82187 7.3±5.5 
2452203.63065 -1.8±18.7 
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Table 4. Radial velocities measured from the 
Na I AA 8183.27,8194.81 doublet in the Clay/LDSS3 
spectra of the two PCEB candidates SDSS J1047+0523 
and SDSSJ1414-0132. 



HJD 




RV [kms- 


J 


HJD 




RV [kms- 


J 


SDSS J1047+0523 




SDSS J1414-0132 




2454237. 


.54059 


214.7±10. 


.4 


2454237. 


.616516 


190.1±11. 


.0 


2454237. 


.64147 


0.0±20. 


.8 


2454237. 


.625451 


192.7±12. 


.1 


2454238. 


.48000 


-127.5±9. 


.5 


2454238. 


.587532 


-78.8±11. 


.2 


2454238. 


.51267 


-107.0±17. 


.2 


2454238. 


.736750 


-lOO.OilO 


.5 


2454238. 


.55507 


-72.5±17, 


.1 


2454239. 


.576845 


4.8±14, 


.9 


2454238. 


.66262 


240. 7±9. 


.8 


2454239. 


.686966 


160.5±7, 


.0 


2454239. 


.54693 


6.6±19, 


.8 


2454239. 


.752785 


209.9±17. 


.2 


2454239. 


.58773 


-119.4±11. 


.9 


2454240. 


.498312 


189.0±11. 


.8 


2454239. 


.61468 


-129.7±8, 


.8 


2454240. 


.569119 


180.2±8, 


.8 


2454239. 


.64270 


-133.0±10. 


.0 


2454240. 


.579801 


181.0±11. 


.0 


2454240. 


.46501 


-31.9±11. 


.6 


2454240. 


.671845 


37.0±7 


.7 


2454240. 


.51077 


109.5±10. 


.5 


2454240. 


.696022 


8.1±10. 


.4 


2454240. 


.59002 


218.7±10. 


.5 


2454240. 


.753426 


-64.1±11, 


.4 


2454240. 


.64421 


162.7±11. 


.7 


2454240. 


.762373 


-62.6±12, 


.0 


2454240. 


.65411 


163.4±9. 


.6 


2454240. 


.790358 


-108.4±7, 


.9 










2454240. 


.801040 


-119.1±13. 


.8 



